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Abstract. The elastic and inelastic collision rates in a gluon gas are calculated. The symmetry 
factor and the phase space integral are discussed in detail. With a symmetry factor and well 
constrained phase space, the same result can be obtained as that of the full phase space without 
the factor. Such an equivalence is illustrated by analytic and numerical calculations for gg — » gg 
and gg — > ggg processes. 



It has been shown in AdS/CFT that the shear viscosity to entropy density ratio rj/s in 
a physical system has a universal minimal bound l/(4-7r) [T]. A recent perturbative QCD 
calculation [2 [3] for a gluon plasma shows that such a limit rj/s ~ 1/Att can be reached by the 
inelastic process gg f-> ggg (23), and the 23 process is 7 times more important than the gg f-> gg 
(22) process. It is in contradiction with the result by Arnold, Moore and Yaffe (AMY) [H [5], 
where the 23 process only provides about a 10% correction to the 22 process. In our previous 
works we have pointed out that such a discrepancy might arise from a multiple counting 

when carrying out collisional integrals of the 23 process over full phase space with a symmetry 
factor. In this note, we will provide a cross check about the above implication by calculating 
collision rates with and without symmetry factors. We will illustrate the equivalence of two ways 
of carrying out collisional integrals, in the constrained phase space with the symmetry factor or 
in full phase space without the symmetry factor. 

For a hot gluon plasma in the vicinity of thermal equilibrium, the collision rates are defined by 
the detailed balance rate in Boltzmann equation [6|. For the 22 process, due to the symmetry 
between the u-channel and i-channel in the matrix element |Mi2->34| 2 , we may choose one 
channel in the the phase space integration and multiply the result by a symmetry factor of 2. 
For example, if we choose the t-channel and limit the phase space to the near forward angle 
scatterings, where q = P3 — Pi, and qr is the transverse component of q with respect to pi in 
the CM frame, the matrix element reads 
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Here a s = g 2 /{Att) is the strong coupling constant, and we have inserted the Debye mass 
rri£) ~ gT to regularize gluon propagator. If we do not limit the phase space in the integration, 




i.e. we take full phase space including both channels, then the factor of 2 is not needed: 
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\M 12 ^u\cm ~ 727r 2 a 2 ^ 2 ; - 2 . (2) 
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Note that the constraint q 2 ~ is removed because both the near forward and backward 
scatterings have small q 2 ^ but only the near forward scatterings have small q 2 [8]. 

With the matrix element, it is straightforward to calculate the elastic collision rate. Using 
the small q approximation, the integral of d^p^cfip^ can be simplified as 
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(2vr) 4 5 4 ( Pl + p 2 - P3 - p 4 )[- ■ • ] « ^ / d V(k~i • q - k 2 • q)[- • • ]. (3) 



The 5 function can be removed by the integral over q, further integral over pi and p 2 gives the 
result R 22 = 2^(3) ~ 4.105a s T for a boson gas with m? D = Aira s T 2 . For a Boltzmann gas 
with m 2 D = 24/7ra s T 2 , the Bose enhancement factor (1 + f®)(l + f®) is absent, and the result 
is i?22 = 3a s T. To check the above analytic results, we can carry out the collisional integral 
numerically using the matrix element in Eq.([2]) in full phase space. The result is shown in Fig(T] 
by the thick solid line, which matches the dashed line of R22 = 3a s T at a s < 0.02 perfectly. The 
difference between the two curves for large a s is due to the fact that the small q approximation 
doesn't hold for large rri£,. 

Another way of calculating the rate is through the cross section [9] with which the elastic 
collision rate in a Boltzmann gas can be written as, 

N g f d 3 pi d 3 p 1 pi ■ p 2 Q Q 
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R22 = ^ fn \\7 \Y„ r /lfe = N g n < (1 - cose 12 )a 22 > . (4) 
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If neglecting the s dependence in 022, it is easy to get R 22 = 3a s T. With the s dependence, 
the average over the pi and P2 gives the thin solid line in Fig{TJ which agrees with the Xu and 
Greiner's result [2] quite well. 

For the 23 process, the exact tree level matrix element |10l \TT\ [T2] is given, e.g., in Eqs. 
(5,6,12) of Ref. [8]. In the limit of one particular configuration in phase space, namely, 



Pl,P2,P3iP4 are hard (their momenta are 0(y / i)), while q(= P2 — Pa = P3 + P5 — Pi) and 
k = ps are soft (momenta are much smaller than y/s), the matrix element becomes 
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|M 12 _345| 2 cm - 3456^— — —j- rf . (5) 

(fcf + m 2 g ) [q^ + m z D ) (k T - q T ) + mf, 

It can be seen that the Gunion-Bertsch (GB) formula |13j is recovered after taking mo,m g — > 0, 
and in this limit, the denominator can be written in a symmetric form Pj^Pj^Pj^. There are 6 
different configurations in the full phase space, represented by the permutation of {3, 4, 5}. If the 
phase space of a single configuration can be isolated and integrated in collisional integrals, the 
result times 6 should be equal to the full phase space result. We will use the above configuration 
where k§ is soft to calculate the rate of the 23 process with the symmetry factor. We now take 
some approximations as follows: (1) For small k, the Bose enhancement diverges as 1 + ~ 1/k. 
(2) As the contribution from small k is enhanced by the divergence, we may assume k <C q and 
|ky| <C |qr|) then the integrals over k and q are decoupled. (3) The recoil from soft gluon to 
hard gluons is negligible, we have 5 4 (pi + p2 — P3 — Pa — Ps) ~ 8 4 (pi + p2 — P3 — Pa), and Eq.Q 
is still valid. So the inelastic collision rate can be written as 
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Note that symmetry factor 6 is added, and the GB formula for the matrix element in literature 
is used. Although complicated enough, most of the integral can be integrated out. We will end 
up with three entangled angular integrals which have to be computed numerically. The final 
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result gives R23 = 4.818a s . Note that the divergence of Bose distribution in soft regime plays 
a very important role. In a Boltzmann gas, these approximations do not hold and the analytical 
approach do not work. In a general case, we should integrate over full phase space without the 
symmetry factor by solving 5 4 (pi +P2 — P3 — Pa — Ps) exactly and numerically. The results from 
two methods are shown as the black solid line and red dashed-dotted line in Figj2j The sizable 
difference at large a s comes from the fact that the small q^ approximation does not hold for 
large mo- 

Another way to regulate the IR singularity is to use the Landau-Pomeranchuk-Migdal 
(LPM) effect |14} [15] to provide a cutoff for the bremsstrahlung gluon. One can simply use 
0(k^ — Ek/lmfp) i n the collisional integral, where the the mean free path can be estimated by 
the inverse of the collision rate l m f p ~ 1/R = l/(i?22 + 2.5i?23). The semi-analytic method is 
still workable and gives the self-consistent equation i?2 3 PM = 6.52a 2 / '(R22 + 2.5i?23 PM )> whose 
solution is R23 PM = + ^5. 2a 2 —R22) ~ 0.99a s T. We can also get a self-consistent result 

by iterating the numerical integration over the full phase space, in which the LPM condition 
applied to the softest gluon for each configuration satisfying the energy-momentum conservation. 
The results from the two methods agree with each other perfectly at small a s as shown by the 
green solid and black long-dashed line in Figj2j 

In summary, we have calculated elastic and inelastic collision rates in a gluon gas with different 
methods to show the significance of the symmetry factor in phase space integrals. The current 




Figure 2. The elastic and inelastic collision rate for a boson gas. Solid lines stand for the 
numerical results in full phase space without the symmetry factor, long-dashed, short-dashed 
and dashed-dotted line stand for semi-analytic results in the constrained phase space with the 
symmetry factor. 



study confirms our previous result about the role of the symmetry factor in the shear and bulk 
viscosity [6j \7\ [8] . We have illustrated the equivalence of carrying out collision integral in the 
constrained phase space with the symmetry factor or in full phase space without the symmetry 
factor. Our result shows that the 22 process is more important than 23 process with reasonable 
coupling constant where the perturbative QCD is reliable. With larger coupling constant, for 
example, in the QGP created at the RHIC/LHC energies, a s ~ 0.2 — 0.3, non-perturbative 
effects should be taken into account to explain the observed small shear viscosity rj/s ~ 1/(4tt). 
The current scheme can be applied to calculate some other thansport coefficients, such as the 
dragging coefficient of heavy quarks, the jet quenching parameter, etc., on which we will look 
at in the future. 
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